COMPLETE GLAUBER CALCULATIONS OF REACTION 
AND INTERACTION CROSS SECTIONS IN LIGHT ION 

COLLISIONS 



I. S. NovikoAfl 

Department of Physics and Astronomy, Western Kentucky University, 
1906 College Heights Blvd, #11077, 

^ BowUng Green, KY 42101-1077, USA 

O 
(N 

Yu. Shabelski 

D 

Petersburg Nuclear Physics Institute, NCR Kurchatov Institute 
^ Catchina, St. Petersburg 188350, Russia 

Abstract 

We calculate the differences between reaction and interaction cross sections in the collisions of 

O 

^ relativistic light ions with A < 40 in the framework of Glauber theory. Although, in the optical 

approximation of Glauber theory these differences are approximately 1% of the reaction cross 
^ sections, they increase up to 3-4% when all scattering diagrams of Glauber theory are included in 

o 

calculation. 
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I. INTRODUCTION 



Various nuclear targets were used to measure interaction cross-sections with stable and 
unstable light isotopes P],E|- The obtained data was used to determine parameters of the 
nuclear density distribution. To do so, the Glauber theory was often used, see [3H5]. The 
reaction cross sections were calculated in the framework of the Glauber theory and compared 
with the experimental data on the interaction cross sections. The values of parameters for 
nuclear matter density can be found in [BHS]. 

The conducted in [U [21 E] analysis relied on the fact that the difference between interaction 
and reaction cross section is negligible. Previously, the difference between o"^^ and cr^ 
has been estimated to be less than a few percent, see flUl ITT] . In [12] authors analyzed 
the reaction and interaction cross section using the "black-sphere" model and rectangular 
uniform distribution of nuclear matter for nuclei with A < 80. It was shown that the 
difference in cross sections is approximately 60 mb, or 4-6% of the reaction cross section 

"ab- 

Since the "black-sphere" model does not proved adequate framework for the calculation 
of reaction and interaction cross sections, we compute the difference between cx^^ and cx^ 
using exact expressions of the Glauber theory as well as expressions obtained in the optical 
approximation. 

II. ELASTIC AND INELASTIC NUCLEUS-NUCLEUS SCATTERING IN THE 
GLAUBER THEORY 

In this section we provide expressions for the reaction, cr^, and interaction, o"^^, cross 
sections and the difference between them obtained in the framework of the Glauber theory. 
Detailed derivation of these expressions can be found in [13]. 

A. Reaction cross section Oj^^ of light ion collisions 

In the Glauber theory, the reaction cross section can be derived using expression for the 
elastic scattering amplitude. As a result we obtain 




(1) 
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where 

SAB{b) = {A\{B\i^ll[l-TN,,{b + u,~Sj)]^ \B)\A) (2) 



with 



TNNib + u,- s,) = ^ [d'q e-^^(^+^-^^)/^^(g), (3) 



fNNil) is the amphtude of elastic nucleon-nucleon scattering and Ui and Sj are the transverse 
coordinates of nucleons. 

In the discussion presented below we used the standard assumptions of the Glauber theory 
(see for example in [H]). First of all, we consider the phase of A-B scattering to be equal 
to the sum of the phases for the scattering of the nucleons of the A-nucleus on the nucleons 
of the B-nucleus. Secondly, we assume that the 3-dimensional distribution of nucleons in 
projectile and target nuclei can be expressed as a product of normalized 3-dimensional single- 
particle densities Ui^z) and p'^i^j, Sj^z), where Ui^z, Sj^z are z-components of the nucleon 
coordinates 

A 

Pj^{ui,Ui^z, ■ ■ ■ ,Ua,Ua,z) = Y\_PA{Ui,Ui^z); 

1=1 
B 

p'BiSl,Si^z---,SB,SB,z) = JJp'i?(Sj,Sj>); (4) 

i=l 

J3„ / 



d-'rip'{r,i) = 1. 

Since we consider forward scattering, we can integrate over longitudinal coordinates Ui,z and 
Sj^z- III this case the expression for Sab as a function of impact parameter b can be written 
in the following form 

»„A r ^ b a b 

SAB{b) = ■■■ Y\pA{ui) <Y]}^~^NN{b + Ui-Sj)]>Y\pB{Sj)Y\dUiY\dSj, (5) 

J J i=i I i,j J j=i i=i j=i 

where pA [ui) and ps (sj) are 2-dimensional single-particle distributions of nucleons in pro- 
jectile and target nuclei. 

In optical approximation, the integral form of Sab can be simplified: 

S%{b) ^ eM-Tor,tm, (6) 

where 



tot r 

ToAb) = ^J dXd%TA{b,)TB{b,)exp 



{b + bi- b 
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(7) 



with 



TA{b) r dzpA (VWT^) . (8) 



B. Interaction cross section of light ion collisions 

The interaction cross sections is defined as a cross section of aU processes which do 
not include target B- nucleus excitation or disintegration (these states are denoted as B*). 
Therefore, the difference between interaction and reaction cross sections can be written as 

<^f^ - = <^iAB ^ AB*) (9) 

Using framework of the Glauber theory, one can obtain following expressions for the 
a{AB AB*): 

(TAB-,AB* = j [JAB{b) - -^^^(fe)]^^^. (10) 

where 

JAB{b) = {A\{B\ I J][l - r^Nib + Ui - sj)] i \A) X 

AA\ jnil - T^^{b-ru, - .;.)]| \B)\A). (11) 
Jab can be simplified using optical approximation 

jZ ^ exp (-r;,(6)) , (12) 

where 

T;,(6) = / d%,d%,TA{b,)T^{b,)e-''^ X 



X ^1 - ] d%TMe ^ j . (13) 

We also provide expression for the cross section of the AB A*B* processes. It can be 
obtained in the same manner as one for a{AB — > AB*): 

CTAB^A*B* = J [lAB{b) - S\j,{b)]d%. (14) 



where 



lAB{b) = {A\{B\ll[[l-rNN{b + Ui-Sj)]\ \B)\A). (15) 



In the optical approximation it reads: 

I%{b) ^ exp (-T;* (6)) , (16) 

where 

T:;Ab) = dXd%,TA{h)Ts{h) i^a'^%e — - 2a^^e —j . (17) 

III. NUMERICAL CALCULATIONS OF cr^^) AND a^^^ DIFFERENCE 

Use of Monte Carlo (MC) integration technique for numerical calculation of [s] was first 
proposed in [151 [16]. MC integration was used for cross section calculations in [T7I - I20] . 

To evaluate Sab, the set of randomly generated nucleon coordinates is required. The 
standard MC approach calls for set of coordinates uniformly distributed in the interaction 
region. That leads to a significant loss of accuracy and an increase in computational time 
due to the fact that often several nucleon coordinates give negligible contribution to the 
Sab- 

The use of Monte Carlo Markov Chain (MCMC) algorithm for generation of nucleon 
coordinates Ui and Sj distributed according to density distributions pA^Ui) and Pb^Sj) proved 
to be a more efficient way to calculate Sab - The Metropolis-Hastings algorithm [21] was used 
to obtain a sequence of random numbers from a defined distributions pA and ps- Generated 
set of coordinates were used to calculate average value of Sab- We apply the same technique 
to (Jab^a*b* and cab-^ab* calculations. 

It is common to use the following parametrization for the amplitude of the elastic nucleon- 
nucleon scattering f^^{q) 

/^^(^) = ^expf-VV (18) 



47r ' V 2 ^ 

where is the total nucleon-nucleon cross section, and {3 is the slope parameter of the 
differential nucleon-nucleon cross-section dependence on q^. We neglect the real part of 
In Nil) since it gives negligible contribution to the reaction cross section. The parameters 



of NN elastic scattering amplitude Eq. ( 18 ) at energy about 1 GeV were taken as 



a^*^ = 43 mb , = 24.3 mb , /3 = 0.35 fm^ (19) 
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The results of the nuclear cross section calculations depend upon the shape of the nuclear 
density distribution, see P]. All calculations presented below were done with Woods-Saxon 
distribution, [22], for projectile and target nuclei: 

PA{r) - 



1 + exp ((r — c)/a) 

In all calculations the value of the parameter a was set to be equal 0.54 fm. 

The c-parameter of Woods-Saxon distribution was varied to match calculated value for 
the reaction cross section o"^, with the experimental value for interaction cross section aixp 
presented in [H [2] . Reaction cross section cr^ was calculated using the expression for Sab 
obtained in optical approximation Eq. [6] and complete expression of the Glauber theory 
(Eq. |2). Experimental cross section, dexp, as well as the values of parameter c calculated in 
optical approximation (OA) and using complete expressions of Glauber theory, are presented 
in Table [T| Obtained results for parameters of Woods-Saxon distribution are in agreement 
with the values for Rrms provided in |9J. 

TABLE I: Experimental interaction cross section fiixp ([1112]) and calculated in optical approxima- 
tion (OA) and in complete Glauber Theory values of Woods-Saxon parameter c (in fm). 



Nucleus (Texp, nib cqa, fm cqi, fm 



Ql2 


853 ±6 


0.85 


1.39 


Ni4 


932 ±9 


1.40 


1.75 


016 


982 ±6 


1.50 


1.90 


pl9 


1043 ± 24 


1.70 


2.10 




1136 ±72 


1.88 


2.25 




1327 ± 14 


2.43 


2.75 



Calculated parameters (cqa and cqi) were used to calculate o"a_b->ab* and aAB-^A'B* 
cross sections for the corresponding stable isotopes. Calculated cross sections <yAB^A*B* 
and (Tab^ab* and computational errors are presented in Table [TTj 

The difference between the values for reaction cross section ct^ab obtained in optical ap- 
proximation and using complete expression of Glauber theory differ by approximately 10%, 
see in pj. Cross sections of the AB — ^ AB* process also depend on the used approximation. 
However, values of cab^ab* obtained by using the expression of the Glauber theory are 
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TABLE II: Calculated in optical approximation (OA) and in complete Glauber Theory aAB~>AB* 
and aAB-¥A* B* cross sections (in mb) 



Nucleus Optical Approximation 


i Ctj Ll I—* V-> A 


Theory 




o AB^A* B* 




u AB^A* B* 


AB^AB* 




169 ±2 


10.4 ± 1.0 


102 ± 1 


34 ±5 




180 ±2 


7.7 ± 1.0 


107 ±2 


34 ±6 


016 


185 ±2 


7.9 ±1.0 


110 ±2 


33 ±4 


pl9 


188 ±3 


6.9 ±1.3 


115 ±3 


36 ±3 




192 ±2 


6.8 ±1.4 


118 ± 1 


41 ±6 


0135 


208 ±4 


5.2 ±1.0 


128 ±2 


42 ±5 



significantly higher than cross sections calculated in optical approximation. 

Since in AB — )■ AB* collisions the projectile nucleus stays in its ground state, we conclude 
that these processes are peripheral. In other words, collisions with an impact parameter 
approximately equal to the sum of the projectile and target nucleus radii contribute more 
to the cab^ab* cross section than collisions with small impact parameters. To illustrate 
this fact, the contribution of different impact parameters to the reaction cross section cr^^, 
aAB^AB* and <7ab^a*b* are shown in Figure 1. All contributions to cross sections were 
calculated for the ^^C-^^C in optical approximation. The error of the MCMC calculation for 
the (Tab^ab* cross section is shown on the graph. 

Obtained results are in agreement with the ones presented in [TU], where contribution 
of impact parameters to and (Tab^ab* cross sections were calculated for ^^0-Al and 
^^0-Pb reactions at 200 GeV per nucleon. 

For the peripheral process we expect the following dependence of aAB^AB* on the atomic 
weight of the projectile nucleus: 

ctab^ab* ~ Ra ~ A^/^, 

when A ^ 1. Therefore, the ratio between cab^ab* and cr^^ is expected to depend on 
atomic weight as 

R* (A) = aAB^AB' /fxS ~ A-- (20) 
with n ~ 1/3. The A-dependence of ratio R*{A) calculated in optical approximation and 
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FIG. 1: Contribution of different impact parameters to the total reaction cross section (solid 
green), to the (Jab-¥A*b* cross section (dashed blue) and to the aAB^AB* cross section (red bars) 
for ^^C-^^C interaction (in semi-log scale). 

using complete expressions of the Glauber theory are presented in Figure 2. 

IV. CONCLUSION 

The differences of cr^ and cr^^ were calculated for light (A = 12-35) stable isotopes 
interacting with C^^ target at approximately 1000 MeV per nucleon. The calculations were 
done in the optical approximation and using complete expressions of the Glauber theory. 

The cross sections (Tab^ab* obtained in optical approximation are approximately 1% of 
the reaction cross section. It is in agreement with the results presented in [11]. However, the 
values calculated using complete expressions of the Glauber theory are approximately 3-4% 
of the reaction cross section. The absolute values of the ctab^ab* calculated using complete 
expressions of the Glauber theory are approximately 30—40 mb. 

When the Glauber theory is used to analyze experimental data on interaction cross sec- 
tion, it is necessary to correct the value of the calculated reaction cross section as 



This leads to the increase of the projectile the -Rrms radii by approximately 2% or 0.05 fm 
for nuclei with atomic weight A < 40. This difference can be additionally increased if the 




'^AB ^ '^ab-^ab*- 



(21) 
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FIG. 2: The values of the R* calculated in optical approximation (blue diamonds) and in complete 
Glauber theory (red circles) as the functions of projectile atomic weight. Data fitting with aA~^ 
function, where n = 1/3, is shown as dashed line. 

projectile nucleus has exited states [2]. 

Since AB — ?■ AB* processes are very peripheral, it is possible that the difference between 
interaction and reaction cross sections increases in the nuclei with halo or skin structures. 
Hence, we believe that when the parameters of nuclear halo and skin are extracted from the 
experimental data, the difference between interaction and reaction cross sections should not 
be neglected. 

This work was supported in part by grant RFBR 11-02-00120-a. 
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